Electron delocalization in bilayer graphene induced by an electric field 



O 
O 

(N 

a 



c3 



c 

o 

> 

o 
o 

(N 
O 
oo 
O 



X 



Mikito Koshino 

Department of Physics, Tokyo Institute of Technology, Meguro-ku, Tokyo 152-8551, Japan 

(Dated: January 16, 2009) 

Electronic localization is numerically studied in disordered bilayer graphene with an electric-field 
induced energy gap. Bilayer graphene is a zero-gap semiconductor, in which an energy gap can be 
opened and controlled by an external electric field perpendicular to the layer plane. We found that, 
in the smooth disorder potential not mixing the states in different valleys (K and K' points), the 
gap opening causes a phase transition at which the electronic localization length diverges. We show 
that this can be interpreted as the integer quantum Hall transition at each single valley, even though 
the magnetic field is absent. 



Since the experimental discovery of monatomic 
graphene [J, 0, 0] , the electronic properties of graphene- 
related materials have been extensively studied. The 
graphene bilayer, which is also experimentally available 
@i HI HI; was shown to have a unique band structure 
distinct from monolayer, where the conduction and va- 
lence bands with quadratic dispersion touch at K and 
K' points in the Brillouin zone [6(. The transport prop- 
erties of the bilayer graphene have been studied experi- 
mentally, \4, M, Lsfl and also extensively studied in theo- 
ries 0, Inl fill [HI . The quantum correction to the 
conductivity, which is important in the low temperature, 
was studied for monolayer graphene [HI, [T3| and for bi- 
layer. 0, E3] A unique property of bilayer graphene is 
that an electric-field applied perpendicularly to the layers 
opens an energy gap between the electron and hole bands 
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20j j . The transport property of the gapped 
bilayer graphene was also studied. 2l| The electric-field 
induced energy gap was observed in recent experiments, 
id 

When the impurity potential is smooth compared to 
the atomic scale and the intervalley scattering is neg- 
ligible, we can treat two valleys (around K or K') as 
independent subsystems. The sub-Hamiltonian within 
each valley generally has no time-reversal symmetry in 
itself, since the time-reversal counterpart of the states 
at K exists at K' . [3, EH In bilayer graphene with 
a gap, interestingly, each sub-Hamiltonian has non-zero 
Hall conductivity even in zero magnetic field. The mono- 
layer graphene can also have an energy gap and the cor- 
responding Hall conductivity when the sub-lattice sym- 
metry is broken [ID, [5^|, while it cannot be controlled 
externally. [23J 

The single-valley Hall conductivity can never be di- 
rectly observed, since it is exactly canceled by the con- 
tribution from the other valley. In this paper, however, 
we find that this quantity strongly influences the elec- 
tron localization properties, and this may be observed in 
the longitudinal conductivity. We show that, when the 
system is under a smooth random potential not mixing 
valleys, the perpendicular electric field causes a phase 
transition at which the localization length diverges. We 
find that this is interpreted as the quantum Hall tran- 
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FIG. 1: Atomic structure of the bilayer graphene. +A and 
—A represents the potential of the top and bottom layers, 
respectively. 



sition at each single valley, where the single-valley Hall 
conductivity changes from one integer to another. To 
demonstrate this, we numerically calculate the localiza- 
tion length in disordered bilayer graphene with electric 
fields, to actually show that the localization length di- 
verges at a certain field amplitude. We then calculate the 
single- valley Hall conductivity and show that the delocal- 
ization is actually associated with the transition between 
different quantum Hall phases. 

Bilayer graphene is composed of a pair of hexagonal 
networks of carbon atoms, which include A and B atoms 
on the bottom layer and A' and B' on the top. As shown 
in Fig, HI both layers are arranged in the Bernal stacking, 
where A atoms are located directly below B' atoms. [f| 
The low energy spectrum is given by the states around 
K and K 1 points in the Brillouin zone. Neighboring A 
and B' sites are coupled to create high-energy bands, 
and remaining A' and B sites form the low-energy bands 
touching at the zero energy. The effective low-energy 
Hamiltonian around -fT-point reads [f| 
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which operates on a wave function (Fj£,, Fjf), where F 
represents the value of the envelope function of K point 
at site X. Here A and — A are the external electrostatic 
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potential at the top and bottom layers, respectively, due 
to the applied electric field, and m* is the effective band 
mass defined by m* = r yi/(2v 2 ), where 71 is the cou- 
pling parameter for the vertical bonds between A and B' 
atoms, and v is the velocity of the monolayer graphenc. 
Q We assume that A and vk are much smaller than 71 
and neglect the terms more than the third order in A/71 
and vk/ji. The band structure of a real bilayer graphene 
is trigonally warped due to the extra coupling parame- 
ter between A' and B, [f3, E| but is neglected here for 
simplicity. The eigen energy of Eq. {!]) is given by 
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k 2 ,. The energy gap extends 



with s = zb and k = w k 2 , 

between e = ±A. The effective Hamiltonian for K' can 
be obtained by exchanging k x + ik y and k x — ik y , giving 
the basically same spectrum. In absence of A, the density 
of states (DOS) becomes constant, p = m*/(2irh 2 ) per 
valley and per spin. 

For the disorder potential, we assume that the length 
scale is much longer than the atomic scale, and neglect 
the intervalley scattering. This should be valid as long 
as the phase coherent time is smaller than the intervalley 
scattering time. [3, EH We also assume that the dis- 
order length scale is shorter than the typical wavelength 
of 27r/fc with k being the wave number from K or K 1 
points, so as to be modeled as a short-ranged potential 
within the valley decoupled Hamiltonian. This is then 
expressed as @ 
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We assume an equal amount of positive and negative 
scatterers Ui — ±u and a total density per unit area n; mp . 
At A = 0, the energy broadening T = H/(2t) becomes 
independent of the energy in the weak disorder limit, and 
is expressed by [§] 
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This will be used as the energy scale characterizing the 
disorder strength. We also introduce unit wave number 
k as H 2 k 2 /(2m*) — T/2 and unit length A = 2ir/k . 

The single-valley Hamiltonian with disordered poten- 
tial TL K = Hq + V belongs to the unitary symmetry class 
when A is non-zero. It becomes the orthogonal class only 
at A = 0, as we have the relation a x 7i K a x — (H K )* with 
the Pauli matrix a x . This is an effective time-reversal 
symmetry within a single valley and should be distin- 
guished from the real time-reversal symmetry connecting 
K and K 1 which always exists at any values of A. 



The single- valley Hamiltonian has a non-zero Hall con- 
ductivity when A / 0. This is estimated by Kubo for- 
mula, 
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and v v are the 



where S is the area of the system, v. 
velocity operators, S is the positive infinitesimal, /(e) is 
the Fermi distribution function, and |a) and e a describe 
the eigenstate and the eigen energy of the system. 

By applying this formula to the ideal i^-point Hamil- 
tonian H.q at zero temperature, we obtain 
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with sgn(x) = The result for K' point is given by 

a xy — —<J xy , so that the net Hall conductivity is exactly 
zero as it should. a xy at the energy gap is quantized into 
different values in the negative and positive A's. When 
Ep is fixed to zero and A is continuously changed from 
negative to positive, we have the Hall conductivity change 
from a xy = —e 2 /h to e 2 /h, and a xy = e 2 /h to -e 2 /h. 

The jump of the Hall conductivity can be intuitively 
explained by the consideration of the energy spectrum in 
the magnetic field. In a uniform field B, the Landau level 
energy at K is obtained by substituting k in Eq. |T]) with 
k + eA/ft with the vector potential A = (0, Bx). [6] This 
results in 



e., n = s^{hw c ) 2 n{n - 1) + A 2 (n > 2) 

£ = £l = -A. (7) 

where oj c — eB/m* and s = ± represents the electron 
and hole bands. The energy spectrum is plotted as a 
function of A in Fig. [21 We have doubly-degenerate Lan- 
dau levels £0 and £1 at E = —A, which moves from the 
electron band down to the hole band as A increases. 
When the Fermi energy is fixed at Ep = and A is 
changed from negative to positive, those two Landau lev- 
els cross Ep, and therefore give a change of a xy by 2e 2 /h, 
twice as large as a Hall conductivity quantum. It is im- 
portant that this discontinuity in a xy is not dissolved in 
the limit of B — > 0, and is actually consistent with the 
discussion above in zero magnetic field. For K' point, we 
have £0 = £1 = +A, leading to the opposite sign of the 
Hall conductivity change. 

Now we numerically calculate the electronic states in 
the disorder potential by exactly diagonalizing the Hamil- 
tonian matrix. We consider a finite square system with 
L x L described by the -fT-point Hamiltonian H K — 
H.q + V, imposing a boundary condition with phase fac- 
tors exp(z0 2; ) and exp(i<f> y ) for x- and y-directions, re- 
spectively. To make the matrix finite, we introduce the 
/c-space cut-off k c = 6/cq, which corresponds to e c = 1ST. 
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FIG. 2: Energy spectrum Eq. (0 plotted against A, at K 
valley in the bilayer graphene with a uniform magnetic field. 

To investigate the localization property, we calculate 
the Thouless number g 7 which is the ratio of the shift 
AE of each energy level due to the change of the bound- 
ary condition, to the level spacing \L 2 p\ with p being 
the density of states per unit area. [2J] We estimate the 
energy shift by AE = ir{\d 2 E{<f) / d(j> 2 \) , where E(<f>) rep- 
resents the eigen energy as a function of the boundary 
phase <j> — <j> x with fixed cj> y , and ( } represents averaging 
over different levels in a small energy region around the 
energy e in question. The localization length L\ oc is es- 
timated by fitting the results to g(L) oc exp(— L/L\ oc ). 
When g •< 1, g is approximately related to the longi- 
tudinal conductivity a xx by a xx ~ (e 2 /h)g. Q We also 
calculate the Hall conductivity by substituting the eigen- 
states of the disordered system to Kubo formula, Eq. (0 . 
For every quantity, we take an average over a number of 
samples with different configuration of the disorder po- 
tential and boundary phase factors 4> x ,4> y . 

The absolute value of a xy generally depends on the 
fc-space cut-off. This is because a xy can be expressed 
as the summation of the contribution from all the oc- 
cupied states below the Fermi energy, unlike <r xx which 
only depends on the states at the Fermi energy. In the 
clean limit, the cut-off at e = ±e c leads to overall shift 
of cr X y(e) by a constant, (e 2 /h)A/e c , which is the contri- 
bution from the missing states out of the cut-off. a xy in 
the disordered system has a similar shift which vanishes 
in the limit of e c — > 0. However, as long as e c is much 
larger than the energy broadening T, it hardly changes 
the function form of a xy {e) around e = except for this 
small constant. 

Figure [3] (a) shows the DOS in disordered bilayer 
graphene, with various A's at fixed T. The energy axis is 
scaled in units of Y. As A increases, the DOS at the zero 
energy monotonically decreases as expected, while the 
disorder potential fills the energy gap with some states. 



(a) Density of states (b) Thouless number (E=0) 
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FIG. 3: (a) Density of states as a function of energy, in bilayer 
graphene with several A's at fixed F. (b) Thouless number g 
at zero energy plotted against A. 
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FIG. 4: Inverse of the localization as a function of energy in 
bilayer graphene with several A's. 

In Fig. [3] (b), we plot the Thouless number g at zero 
energy as a function of A, for several system sizes. We 
naively expect that the conductivity at zero energy mono- 
tonically decreases in increasing A since the DOS be- 
comes smaller, but this is not the case here. At every 
system size, g goes up first when A starts from 0, and 
then changes to decrease at A ~ 0.8r. Significantly g be- 
comes independent of the system size around the peak, 
while in other region it drops exponentially as the system 
size increases. This suggests that the electronic state at 
zero energy is localized both in small and large A's, but 
delocalized only at A ~ 0.8r. We calculate g also for 
other energies and obtain the localization length L\ oc as 
a function of energy from its size dependence. In Fig. |4j 
we plot \/L\ oc against energy for several values of A. As 
A increases from zero, l/L\ oc gradually decreases, and 
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FIG. 5: Hall conductivity of a single- valley (K) Hamilto- 
nian of bilayer graphene with several A's, plotted against the 
energy. 
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FIG. 6: Hall plateau diagram against A and energy E. 
Dots represent the critical energies estimated from the size- 
dependence of a^ y . Dashed curves are guide for eyes for the 
phase boundary. 



the value around s = reaches almost zero (i.e., the 
states extended) at A ~ 0.8I\ The center value recovers 
to non-zero again in even larger A, and the peak grows 
rapidly as A increases; the states at the gap region be- 
come localized more and more. 

The corresponding plots for the single- valley Hall con- 
ductivity a^y are shown in Fig. [5] a^ y becomes com- 
pletely at A = (not shown), and the value generally 
goes up as A increases. At each single A, we notice that 
a^y moves in a specific direction as the system size L ex- 
pands. At A = 0.4r, cr£y decreases in all the energy as 
L increases. At A = 1.6r, on the contrary, a^ y around 
e = evolves in the upward direction, oppositely. The 

critical value of A where the a^f„ switches its evolving 

* y 



direction depends on the energy. At E = 0, this is close 
to A/r = 0.8. 

In the integer quantum Hall effect, the Hall conduc- 
tivity in the infinite system is quantized as long as the 
states at the Fermi energy are localized. [27[ In finite- 
size systems, the Hall conductivity averaged over disorder 
configurations generally takes non-integer values, while it 
gradually approaches the quantized value as the system 
size expands. |28j Therefore the quantized Hall conduc- 
tivity in the infinite system can be estimated from the 
scaling behavior of finite-size values. In above results, we 
expect that the region of decreasing a^ y becomes the Hall 
plateau with <j^ y = 0, and that of increasing <j^ y becomes 
the plateau with e 2 //i, in an infinite system size. The 
point where the a^ y switches the scaling dependence can 
be regarded as the phase transition between the different 
Hall phases a^ y = and 1. Significantly this point coin- 
cides well with the divergence of the localization length 
L\ oc in Fig. [4j The feature is most striking at A = 0.8r, 
where a^ y at e = is about to change its size dependence, 
and L\ oc at the center is also near divergence. This is a 
consistent result since the electronic states must be de- 
localized at the critical point separating different Hall 
plateaus. 
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Indeed, when the energy gap width is much larger than 
the disorder strength, i.e., A/r 3> 1, we reasonably ex- 
pect that a^ y around zero energy is quantized at e 2 /h, the 
mid-gap value in the clean-limit, since the gap remains 
almost intact in the weak disorder. When the energy gap 
is small enough that A/r <C 1, on the other hand, it is 
natural that a^ y vanishes to zero as in the zero-gap case, 
since the effect of A is completely washed out by the 
strong disorder. It is then inevitable to have the phase 
transition at the intermediate value of A/r, which is pre- 
sumably of the order of 1, while the factor may depend 
on the specific disorder model. 

At A > 0.8, we have two critical energies which sep- 
arates upward and downward moving region in a^ y in 
Fig. [5l The localization length in Fig. [4] is huge around 
the corresponding energies, so that it becomes harder to 
specify the critical points out of them due to the numer- 
ical error. Fig. [6] shows the phase diagram speculated in 
the infinite system, where each phase corresponds to the 
Hall plateau with quantized a^ y . We determine the phase 
boundary by taking points where a^ y changes its moving 
direction. The dashed curves separating the phases are 
guides for eyes. The phase diagram is symmetric with 
respect to A = 0, while the sign of the Hall conductivity 
is opposite between negative and positive A's. 

We cannot eliminate the delocalized states as long as 
K and K' valleys are decoupled, in the following sense: 
In the clean system, as previously discussed, the single- 
valley Hall conductivity at the gap is quantized at 
— 1 and 1 (in units of e 2 /h) at negative and positive A, 
respectively. In the disordered system, we expect that 



5 



<7jr at the gap remains quantized when the gap width is 
large enough compared to the energy broadening due to 
the disorder. Thus the Hall conductivity at a fixed Fermi 
energy definitely changes from —1 to 1 when A changes 
from — Aq to +Ao with sufficiently large Ao. This re- 
quires that the delocalized states appear at that energy 
somewhere between — Ao and +Ao, because otherwise 
the Hall conductivity stays constant. 



27j The present 



calculation indeed suggests that there are two delocal- 
ized points of A (negative and positive) at each fixed 
energy. 

The delocalized states would disappear in presence 
of the short-ranged disorder with atomic length scale, 
which induces intervalley coupling between K and K 1 
and makes the original time-reversal symmetry effective. 
The localization length at the phase boundary is expected 
to diverge as the inter-valley scattering rate goes to zero. 
The trigonal warping effect due to the extra hopping pa- 
rameter, which is also neglected here, does not kill the 
delocalized states since it does not mix the valleys, while 
it g enerally weaken the electron localization at A = 0. 

m 



In this paper we proposed that the external electric 
field applied to the bilayer graphene causes a phase tran- 
sition at which the localization length diverges. The tran- 
sition can be interpreted as an analog of quantum Hall 
transition at each single valley. While the actual Hall 
conductivity is never directly observed, it may be possi- 
ble to observe the divergence of the localization length in 
a temperature (T) dependence of the conductivity a xx , 
as in the conventional integer quantum Hall effect. [29( 
Figure EJb) can be approximately viewed as plots of a xx 
at different T's, when we use a relation a xx ~ (e 2 /h)g, 
and regard L as the phase coherent length at T. When 
the system goes through the phase transition, a xx ex- 
hibits a peak at the phase boundary and its peak width 
becomes narrower as the temperature is lowered. [29| To 
observe this, the temperature must be low enough that 
the phase coherent length exceeds the typical localization 
length lg. Our calculation is based on the a typical dis- 
order model used in theories, while it may not describe 
some specific situations in the real bilayer graphene cur- 
rently available, which is supposed to include puddles 
(long-range inhomogeneity) or impurity-bound states de- 
scribed by the variable-range hopping model. We leave to 
the future work the study on those effects on the present 
phenomenon. 
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